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A theoretical study of the intense-field multiphoton ionization of hydrogenlike systems is per- 
formed by solving the time-dependent Dirac equation within the dipole approximation. It is shown 
that the velocity-gauge results agree to the ones in length gauge only if the negative-energy states 
are included in the time propagation. On the other hand, for the considered laser parameters, 
no significant difference is found in length gauge if the negative-energy states are included or not. 
Within the adopted dipole approximation the main relativistic effect is the shift of the ionization 
potential. A simple scaling procedure is proposed to account for this effect. 



PACS numbers: 32.80.Fb, 31.30.J- 



INTRODUCTION 



Future experiments using an electron-beam ion trap 
(EBIT) at the Linear Coherent Light Source (LCLS) at 
Stanford and X-ray Free Electron Laser (XFEL) at Ham- 
burg are expected to permit the study of photoabsorp- 
tion processes of highly charged atomic ions in the wave- 
length range of 0.1 to 100 nm with a peak intensity up 
to lO^^W/cm^ or even higher. Also at the GSI (Darm- 
stadt) experimental investigations of highly charged ions 
exposed to intense laser fields are planned within the 
SPARC project. The analysis of these experiments will 
require a relativistic treatment of the ion-laser interac- 
tion. 

Clearly, a full treatment demands to solve the time- 
dependent Dirac equation (TDDE) incorporating also 
the spatial dependence of the vector potential. Since 
such a treatment is very demanding, most earlier treat- 
ments adopted simplifications. Low-dimensional mod- 
els are especially popular. Starting first with a one- 
dimensional treatment [l[ , elaborate two-dimensional cal- 
culations have been reported, e.g., in j2-4]. However, 
such models can certainly not provide quantitative pre- 
dictions and it is, at least a priori, not even clear whether 
they are always qualitatively correct. Similarly to the 
non-relativistic case, also simplified ionization models 
like the strong- field approximation 0-11], semiclassical 
tunneling theory [§], or classical models like in (Tol - [l^ 
have been proposed. However, in order to allow for 
quantitative predictions or the validation of such sim- 
plified models a full-dimensional solution of the TDDE 
is needed. 

Very recently, a three-dimensional solution of the 
TDDE for hydrogenlike systems has been reported in 
[is'l. The radial solutions were expanded on a grid and 
the TDDE was either solved by a direct propagation on 
the grid or using a spectral expansion in field-free eigen- 
states. The spatial dependence of the vector potential 
of the carrier part of the laser pulse was also consid- 
ered, while the one in the envelope was ignored. The 
velocity gauge was used and the importance of includ- 



ing the negative-energy (NE) states was emphasized in 
the case of a treatment beyond the dipole approxima- 
tion, even for the considered laser parameters where the 
photon energy is insufficient to produce real positron- 
electron pairs. Based on general theoretical considera- 
tions it was conjectured that in length gauge the impor- 
tance of the NE states may be reduced. On the other 
hand, it was concluded on the basis of the numerical re- 
sults that within the dipole approximation the inclusion 
of NE states is not needed, even if the TDDE is solved 
in velocity gauge. A comparison of the numerically ob- 
tained ionization rates for various nuclear charges with 
the ones in non-relativistic approximation showed that, 
expectedly, increasing differences are found with increas- 
ing charge. However, it was concluded that the ionization 
rate (shown as a function of the peak value of the laser 
field) obtained within the relativistic TDDE calculation 
may be larger or smaller than the non-relativistic result. 
The authors found the higher rate easier to understand 
and could only speculate on possible reasons for the lower 
one. 

A solution of the TDDE within the length gauge was 
reported more recently in il4]. A direct time propagation 
on a grid was used. As in [ij] the electron-nucleus in- 
teraction is described by the unmodified non-relativistic 
Coulomb interaction. While solely the dipole approxima- 
tion is adopted, not only results for one-electron, but also 
for two-electron ions are reported. In the latter case the 
electron-electron interaction is also described by the non- 
relativistic Coulomb interaction. No explicit discussion 
of the inclusion or omission of the NE states is given. 
Comparing for Ne^+, Ne«+, \]'^'^+ , and U^°+ the pho- 
toionization cross sections obtained by solving the TDDE 
for 10-cycle pulses (using a single pulse with one set of 
laser parameters for every ion) with the ones of relativis- 
tic perturbation theory, a good quantitative agreement 
is found. This result is, of course, expected, if the laser 
parameters are chosen in a way that perturbation theory 
is applicable, as was the case in [l3] where relatively low 
intensities (in relation to the ionization potentials and 
photon frequencies) were considered. 

In this work a further approach for solving the TDDE 
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is presented. It is based on the spectral expansion in 
field-free eigenstates where the radial wavefunctions are 
expressed in a i?-spline basis. It may be seen as an 
extension of the corresponding approaches for solving 
the non-relativistic time-dependent Schrodinger equation 
(TDSE) for one- and two-electron atoms in, e.g., 
U^ . for molecular hydrogen [l9l - l23l |. and for in princi- 
ple arbitrary molecules within the single-active-electron 
approximation [23, While and considered 
mostly (or solely) one-photon ionization, an extension to 
multiphoton ionization is presented. Clearly, such cal- 
culations are very demanding as they require larger ex- 
pansions, because more angular momenta are involved. 
Within the dipole approximation a systematic investiga- 
tion of the importance of NE states as well as the con- 
vergence behavior within either length or velocity gauge 
is performed. A simple scaling relation is proposed that 
allows to relate relativistic TDDE solutions to the ones 
obtained with the non-relativistic TDSE, if both are per- 
formed within the dipole approximation. 

In this work atomic units (e = mg = h ~ 1) are used 
unless specified otherwise. 



II. THEORY 

The dynamics of a highly charged atomic ion exposed 
to an external electromagnetic field is considered by 
means of both a relativistic and a non-relativistic treat- 
ment within the dipole approximation. The vector poten- 
tial A(t) is chosen in the form of an iV-cycle cos^-shaped 
laser pulse that is linear-polarized along the z axis. 
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where oj is the radiation frequency and T ~ 2t:N/uj is 
the pulse duration. Ignoring finite nuclear size effects, 
the interaction of an electron with the nucleus may be 
approximated by the Coulomb potential 



U{r) 



(2) 



where Z is the nuclear charge. In all calculations pre- 
sented in this work the hydrogenlike system is initially 
prepared in its ground state. 



Relativistic calculations 



and the interaction with the electromagnetic field V{t) 
can be presented within the dipole approximation either 
in the velocity (V) or length (L) gauge. 



V^{t) ^ca- A{t) = ca;,A{t), 
V^{t) ^r-F{t) ^zF{t). 



(5) 



Here, f3 and the components of the vector a are the Dirac 
matrices, p = — «V is the momentum operator, c ~ 137 
is the speed of light, and F{t) = —dA{t)/dt is the electric 
field. 

As is well known (see, e. g., [1^), the eigenstates of 
can be presented as four component spinors 

where XK.,m(j) is an Is coupled spherical spinor, k is the 
relativistic quantum number of angular momentum, re- 
lated to the orbital and total angular momenta, I and j, 
as 

^ f-(j + l/2) = -(/ + l), forj=/ + l/2, 
" \ J + 1/2 = I, for j=l- 1/2; ^ ' 

and the radial functions P^ir) and Q^if) are solutions of 
the coupled equations 




E]Q^{r) + c 
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r dr 
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r dr 



Qn{r) = 
P.{r) = 0. 



(8) 



For two spinor states, $i and character- 
ized by the quantum numbers {ui, Ki, ji,li,mi} and 
{nf,Kf,jf,lf,mf}, respectively, the time-dependent ma- 
trix element Vfi{t) ~ f\V{t)\^i) can be written as 



Vf^{t) = 6rnf,^. S\lf-h\,l Wf^ Mf,{t) (9) 



where 



Wf. = (-l)«-"/(-iy>+V2^(2^-^ + i)(2j, + 1) 
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jf 1 ji 
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(10) 



and 



M]^,it)^Fit) drr P.,(r)Pf{r) + Q.,{r)Qf{r) , (11) 



The relativistic dynamics of the quantum system is 
governed by the TDDE 



where is the field- free Dirac Hamiltonian i/^, 



ca • p -F + C/(r), 



(3) 



(4) 



Mlit) 



with A 



icA[t) / dr [(l + A/0^*(r)Q/(r) 
-{\-^f;)Q,{r)Pf{r) 



(12) 



In order to describe both bound and continuum states, 
the atom is confined within a spherical box boundary of 
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radius R. This leads to the discretization of the contin- 
uum spectrum, whereas the bound states remain unmod- 
ified, if R is chosen sufficiently large and not too highly 
excited Rydberg states are considered. 

Introducing in a region [0, R] a _B-spline set consisting 
of rt-|-2 spline functions Bi{r) of order k, the radial func- 
tions P{r) and Q{r) can be expanded in a _B-spline basis 
as 



P{r) = J2P^B,+,{r) Q(r) = ^g,S,+i(r) (13) 

1=1 



1=1 



where the first and the last spline are removed from the 
expansion to ensure the boundary conditions P(0) = 
Q(0) = P{R) = Q[R) = 0. Defining the coefficient vec- 
tor 



(Pl,gi,P2,'72,---,P«,gn), 



(14) 



the system of equations ([5]) is transformed into a,2nx2n 
generalized banded eigenvalue problem which is effi- 
ciently solved using LAPACK routine DSBGVX. This 
procedure yields for every value of k exactly n negative 
energy solutions and n positive energy solutions. 

As has been discussed in literature since the first solu- 
tion of the time-independent (stationary) Dirac equation 
using B splines [13], there is the problem of the occur- 
rence of spurious states [1^ that are non-physical. Dif- 
ferent procedures were proposed to avoid this problem 
[29I - I31I . While there can be a problem with their iden- 
tification in the case of many-electron systems, it is the 
lowest positive energy state for k > that represents a 
spurious state in the present case. 



B. Non-relativistic calculations 



The non-relativistic TDSE 



i^ = [H, + V{m{t) 



(15) 



where Hq is the non-relativistic field-free Hamiltonian 

Ho, 



(16) 



and in contrast to the relativistic case, the interaction 
with the electromagnetic field V(t) is given in the velocity 
gauge by 



(17) 



Similar to the relativistic case, the eigenstates of Hq, 
^';m(r) = r~^i?/(r)Yim(f), are obtained by projecting the 
radial function Ri{r) onto the i?-spline basis 



(18) 



and transforming the radial Schrodinger equation into 
the n X n generalized banded eigenvalue problem with 
respect to a coefficient vector 



C = (pi,P2, • • • ,Pn), 



(19) 



which yields n solutions for every orbital quantum num- 
ber I. 

For two states, $i and characterized by the 

quantum numbers {ni,li,mi} and {nf,lf,mf}, respec- 
tively, the time-dependent matrix element Vfi{t) = 

f\V{t)\^i) can also be written in form of Eq. |9l but 
with 



Wf. = (-l)'/-"V(_i)'/y^(2;/ + i)(2/, + l) 



If 1 h 



-nif rrii J \0 



1 /, 



and 



MfM^F{t) / drrR,{r)Rfir), 



Ml{t) = -zA{t) / drRf{r) 



-R,{r) + R[{r) 



(20) 



(21) 



(22) 



C. Time propagation 

Within the spectral approach, the integration of the 
TDDE © and the TDSE ^ is performed by expanding 
the function ^'(i) describing the dynamics of the system 
in the basis of field- free eigenstates ^k, 



*(r,t) = ^C 



.K(<)e-'^-*$^(r) 



(23) 



K 



where the compound index K represents the full set of 
quantum numbers and the coefficients Cxit = —T/2) 
are set to zero for all states except for the initial state for 
which its value is set to 1. Since the quantum number 
TO is conserved for both cases and the ground state is 
chosen as initial state, the value of m is fixed to 1/2 for 
the relativistic case and to for the non-relativistic case. 

Substituting Eq. ^ into Eq. © or TDSE ^ the 
latter ones are reduced to a system of coupled first-order 
ordinary differential equations. 



C'Kit) - 5]e'(^"-^-')* Vkk'(<)Ck'(<) 



(24) 



K' 



which is integrated numerically using a variable-order, 
variable-step Adams solver. 

The ionization yield (after the pulse) is then defined as 



yio„ = ^|CK(t = T/2)|' 



(25) 



K 



where the summation in Eq. ()25p is performed over all 
discretized continuum states. The convergence can be 
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controlled by varying the value of Zmax (or jmax) which 
limits the number of different symmetries involved in the 
summation ([23]). 

In the case of the TDDE the question of a proper treat- 
ment of the spurious states mentioned in Sec lII Al arises. 
The success of the spectral approach relies on the com- 
pleteness of the states included in the summation, at least 
for a given box size. In fact, already in the case of the 
non-relativistic Schrodinger equation the use of a finite 
box leads to the occurrence of non-physical pseudo states 
(see, e.g., In this case it turned out to be in fact 

important to include those states in the spectral expan- 
sion, since otherwise some relevant part of the Hilbert 
space is omitted. In order to decide on the proper treat- 
ment of the spurious states of the Dirac equation (due to 
the _B-spline expansion) , a careful check of the relativistic 
sum rules 128] was performed that also served as a check 
of the proper implementation of the length and velocity 
forms of the dipole-operator matrix elements. On this 
basis it was concluded that the spurious states could be 
omitted from the basis for solving the TDDE, since the 
sum rules were fulfilled in this case. The agreement of 
the TDDE solutions obtained within the length and the 
velocity gauge discussed below is another indication that 
the omission of the spurious states is appropriate. 
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FIG. 1: Convergence study with respect to the maximal value 
of the total angular momentum jmax. The relativistic calcula- 
tions of the ionization yield for an ion with the nuclear charge 
Z — 50 exposed to a 2Q-cycle cos^-shaped laser pulse with 
peak intensity / = 5 x 10^'^W/cm'^ and a photon energy of 
500 a.u. are performed using the length (L) or velocity (V) 
gauge either including or excluding the negative energy (NE) 
states in the expansion (|23|l . 



D. Computational details and scaling of the TDSE 
with Z 



feet, cf. [13|. 



For the sake of consistency, the same B-spline basis set 
is used for both the relativistic and the non-relativistic 
treatment. Typically, the values A; = 9 and n = 500 are 
used to construct an almost linear knot sequence in which 
the first 40 intervals increase with a geometric progres- 
sion using g = 1.05 and all following intervals have the 
length of the 40th one. Such a choice ensures an accu- 
rate numerical description in the vicinity of the nucleus 
and, at the same time, a sufficient completeness for a de- 
scription of the discretized continuum [sl]- Depending 
on the nuclear charge Z, a box size of i? = (250/Z) uq is 
adopted. This choice of R refiects the well-known scaling 
property of the time-independent Schrodinger equation 
of hydrogenlike systems. With the substitution r' = r/Z 
and E' = Z'^E for the position r and the energy i?, re- 
spectively, the eigensolutions for a hydrogenlike system 
with charge Z reduces to the one of the hydrogen atom 
with Z = \. 

Within the dipole approximation this property is pre- 
served also for the TDSE, if the pulse parameters are 
also properly scaled, e.g., the time as t' = t/Z"^, the 
laser frequency as to' = Z'^uj, the amplitude of the vector 
potential as A'q = ZAq, and the laser peak intensity as 
/' = Z^I. Since this property does not persist in case of 
the relativistic treatment, and the properly scaled solu- 
tions of the TDDE are thus not identical anymore, any 
deviation from the non-relativistic prediction based on 
the scaling relations can be classified as a relativistic ef- 



III. RESULTS AND DISCUSSION 

A. Convergence behavior 

The handling of the NE solutions of the Dirac equation 
is an important issue and steered considerable attention 
in the literature, see [ij, [s^ and references therein. The 
question is whether the NE states should be removed 
from the basis as long as the creation of real positron- 
electron pair is energetically out of reach. Such a sit- 
uation occurs, for example, in case of ionization of an 
ion with the nuclear charge Z = 50 exposed to a pulse 
with photon energy 500 a.u. (see Fig. [1]). Whereas only 
3 photons are sufficient for ionization, more than 70 pho- 
tons are required for positron-electron pair creation. As 
is demonstrated in Fig. [U the answer to the question of 
the proper treatment of the NE states depends already 
within the dipole approximation adopted here clearly on 
the gauge and thus the interaction operator used in the 
time propagation. For the length gauge the exclusion of 
NE states has virtually no effect on the final result. In 
contrast, for the velocity gauge the results are obviously 
different. If the NE states are included in the time propa- 
gation, the converged result agrees with the one obtained 
using the length gauge, whereas without the NE states 
the converged result differs in this example by a factor of 
about 1.5. Therefore, these states are absolutely neces- 
sary for obtaining the correct result. 
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From the practical point of view the finding is also in- 
teresting, since the omission of the NE states and thus 
half of the total number of states reduces the numerical 
efforts tremendously. Clearly, doubling the number of 
states leads to an increase of the number of operations 
per time step by a factor 4, since the number of transi- 
tion dipole matrix elements increases by this factor. In 
the present example (and for a given number of jmax) 
the time propagation in length gauge without the NE 
states is by about a factor of 6 faster compared to the 
one where the NE states are included. The additional 
time difference (factor of about 1.5) arises from an in- 
creased number of time steps required for convergence in 
the used adaptive time propagation, if the NE states are 
included. In fact, the length-gauge time propagation by 
itself is found to be about 6 times faster than the one 
performed in velocity gauge, even if the NE states are 
included in both of them. This factor is due to the finer 
time grid needed for convergence in velocity gauge. Con- 
sidering both effects together, even a speed-up by a factor 
of 50 is found when comparing the length-gauge calcu- 
lation without NE states and the velocity-gauge variant 
with NE states that for sufficiently large value of jmax 
both yield practically identical results. 

However, Fig. [T] reveals also that the calculations in 
velocity gauge converge faster with respect to the quan- 
tum numbers j included in the calculation. Therefore, 
the efficiency gain of the length gauge is smaller than 
the value given above. In the concrete example shown in 
Fig. [T]the (within better than 0.1%) converged length- 
gauge calculation (jmax — 15/2, without NE states) is by 
a factor of about 40 faster than the (within about 10~^ %) 
converged velocity- gauge result (jmax = 11/2, with NE 
states). Even compared to the only about 0.5% con- 
verged velocity-gauge result with jmax = 9/2 (with NE 
states) there is still a factor of about 30 in time gain. On 
the other hand, the question of the most efficient choice 
of the gauge depends on the laser parameters. If few- 
photon processes, in the most extreme case one-photon 
ionization at low intensities, are considered, the average 
angular momentum j transferred to the ion is small. In 
such a case the need for a larger value of jmax m the 
length-gauge calculation could even overcompensate the 
gain from the exclusion of the NE states. However, if 
the average number of absorbed photons increases (due 
to a smaller ratio of the photon frequency with respect 
to the ionization potential or due to a higher laser inten- 
sity), the relative increase in j due to the slower conver- 
gence leads to a smaller increase of the total number of 
states than the inclusion of the NE states. In this mul- 
tiphoton regime calculations in length gauge appear to 
be much more efficient. This finding for the TDDE solu- 
tion seems to differs from the one for the non-relativistic 
TDSE where the velocity gauge is often supposed to con- 
verge faster than the length gauge. 

For extremely high intensities around and above the 
critical field strength Fcr = m 2.57 x 10^ where real 
pair production is possible, the inclusion of the NE states 
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FIG. 2: Ratio of the ionization yields that are obtained by 
either excluding or including the negative energy states in 
the expansion (|23p adopting either length (L) or velocity (V) 
gauge. The calculations are performed for an ion with two dif- 
ferent values of nuclear charge Z by varying the wavelength 
A of a 20-cycle cos^-shaped laser pulse with the peak inten- 
sity Z^ X 10^"^W/cm^. The character of the ionization pro- 
cess changes from single-photon ionization (XZ^ = 40 nm) to 
three-photon ionization (XZ^ — 260 nm). 



is, of course, required also in length gauge. In fact, al- 
ready for the parameters discussed in the context of Fig. [T] 
we find a relative deviation between otherwise converged 
length-gauge results with and without NE states of about 
2.2 X 10^^. This appears reasonable, since it is of the 
order of the (reciprocal) rest energy « 5.3 x 10~^. 
For very high laser intensities the gain from omitting the 
NE states is thus lost and future calculations will have 
to show whether in that regime length-gauge calculations 
can still profit from the need of a sparser time grid, or 
whether the faster convergence with j persists in velocity 
gauge and may make this gauge more efficient. 

In |.13i] it was concluded on the basis of the numerical 
results that the inclusion of the NE states is not cru- 
cial for solving the TDDE, if the dipole approximation 
is adopted. This result was found despite the fact that 
the calculation was performed in velocity gauge. There- 
fore, this finding appears to contradict our conclusions. 
However, it turns out that the importance of the inclu- 
sion of the NE states depends in fact on the character 
of the multiphoton ionization process, i. e. on the num- 
ber of photons. Figure [5] shows the ratio of the ioniza- 
tion yields obtained either with or without inclusion of 
the NE states for both gauges and two values of the nu- 
clear charge Z. In agreement with the findings discussed 
above, the inclusion of the NE states has practically no 
influence for the length-gauge results, independently of 
the number of photons involved or the nuclear charge. 
This changes, however, if the velocity gauge is used in 
the time propagation. The importance of the NE states 
increases with the order of the multiphoton process (the 
number of photons required for ionization) and with the 
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nuclear charge. Within the one-photon regime there is 
practically no influence of the NE states. Thus the find- 
ing in [131 is confirmed that even in velocity gauge an 
inclusion of the NE states is not required for calculations 
in the dipole approximation. However, it is to be under- 
stood that those findings apply only to the case where 
one photon is sufficient to reach into the ionization con- 
tinuum. Quite remarkable is also the strong dependence 
of the importance of the NE states on Z. For exam- 
ple, for the wavelength A = 260/Z^nm the exclusion of 
the NE states leads to a three times smaller value of the 
ionization yield for Z = 80, whereas for Z = 40 the de- 
crease is only about 10%. The Z dependence can also 
explain the negligible effect of NE states within the dipole 
approximation found in Fig. 1 of [l3l | where Z — \ was 
used. Furthermore, only the probability of remaining in 
the initial ground state was considered and this quantity 
is found to converge much faster than, e.g., the ioniza- 
tion yield. This could possibly make it also less sensitive 
to the NE states. 




15 20 25 30 35 40 45 50 
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eV 

FIG. 3: Ratios of the relativistic ionization yields to the non- 
relativistic ones (shown in the insert) obtained for 5 different 
values of the nuclear charge Z by varying the carrier frequency 
cj of a 20-cycle cos^ -shaped laser pulse with the peak intensity 
of Z^ X 10"W/cm^ 



B. Single-photon ionization 

As has been discussed in Sec. IIIDi the solution of the 
TDSE for different nuclear charges Z gives identical ion- 
ization yields, if the laser pulse parameters are scaled 
properly. For example, the insert in Fig. [3] shows the 
non-relativistic ionization yields obtained for an ion with 
the nuclear charge Z exposed to a 20-cycle cos^-shaped 
laser pulse with a peak intensity of x lO^^W/cm^ for 
photon energies varying in the range between 15 Z"^ eV 
and 45 eV. Since the ionization potential of the ion 
is equal to 13.6 eV, the ionization should occur via 
absorption of a single photon. In order to study the rel- 
ativistic effects in this one-photon ionization regime, the 
TDDE is solved for the same system and for 5 different 
values of Z (in between 40 and 80 with a stepsize of 10). 

The ratio of the relativistic to the non-relativistic ion- 
ization yield is shown in Fig. [3] Two features can be 
observed. First, above ISZ^eV the ratio decreases with 
increasing photon energy. Thus relativistic effects are 
more pronounced for higher photon frequencies cj. As 
could be expected, this effect becomes more and more im- 
portant as Z increases. Second, a discontinuity develops 
with increasing Z at a photon energy of about 18 Z^ eV. 
Especially for large values of Z the ratio for the lowest 
considered photon energy drops down substantially for 
increasing Z . This is a manifestation of the fact that the 
ionization potential of the ion increases as a consequence 
of the relativistic velocity of an electron in the vicinity 
of strong Coulomb potentials. The shift of the ionization 
potential 

A/p = lf{Z)-ll\Z) = (l - (2^) 

rapidly increases with increasing Z . (For example, A/p 
Z4/(8c2) for smaU Z) Thus, for Z = 80 the absorption 



of a single photon with the energy 15 eV is not suf- 
ficient for ionization anymore and the character of the 
ionization process changes from single-photon to two- 
photon ionization. For the considered intensities the lat- 
ter process possesses of course a much smaller proba- 
bility. In fact, the ratio would be even smaller, if the 
ratio at 15 eV would correspond to pure two-photon 
ionization. The finite spectral width of the adopted 20- 
cycle laser pulse allows, however, one-photon ionization 
to occur even at this energy and increases thus the ion- 
ization yield. Therefore, a finer photon-frequency grid 
in between 15 Z^ eV and 20 Z^ eV would show a much 
smoother behavior than the one visible from Fig. [31 

In [l^l (Fig. 5) relativistic effects were considered by 
a comparison of the ionization rates obtained with the 
TDSE and the TDDE where the latter was solved for 8 
different values of Z. However, there the behavior was 
studied as a function of the laser field amplitude (also 
scaled by the corresponding non-relativistic scaling re- 
lations) and for a fixed photon frequency. The authors 
discuss stabilization, since the ionization rate does not in- 
crease monotonously with the field strength, but instead 
decreases for intensities beyond about Fq = IZ'^. The 
for large values of Z increasing value of Fq for which the 
ionization rate has its maximum is then explained by the 
stabilization criterion of Gavrila (see [sS] and references 
therein) and the conjecture that due to the lower energy 
of the ground state in the relativistic case the condition 
for the occurrence of stabilization shifts to higher field 
strengths. 

In fact, the dependence of the ionization rates on Z 
discussed in Ts'] may be quantitatively understood from 
the scaling relations together with the lowering of the 
ground-state energy due to relativistic effects. This is 
illustrated in the following way. From the condition 
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FIG. 4: Ionization rate (in units of a.u.) of a hydrogenlike 
ion exposed to a monochromatic laser field with the frequency 
tj = as a function of the amplitude of the electric field 
(given in units of a.u.) The black solid curve presents the 
(non-relativistic) Floquet ionization rate, whereas the other 
curves present the ionization rates obtained using the scaling 
relation [28] in order to estimate the relativistic ionization rates 
for three difi'erent values of Z. 



Ip°^{Z) = Ip'-'{Z') one finds a scaled nuclear charge 



Z' 



y2c2(l- v/T^ZV?), (27) 



for a given true charge Z. This allows to estimate the 
relativistic ionization rate from a non-relativistic calcu- 
lation using 



r Fo 



(28) 



if the dependence of the non-relativistic ionization rate on 
the photon frequency is ignored. In order to obtain non- 
relativistic ionization rates r(Fo), Floquet calculations 
were performed as a function of the peak amplitude Fq 
of the electric field using program STRFLO [35|. They 
are shown in Fig. |4]and are in reasonable agreement with 
the TDSE rates in [13]. Based on the non-relativistic Flo- 
quet rates the scaling relation allows to estimate the 
relativistic rates that are also shown in Fig. U for Z = 36, 
54, and 86. The qualitative agreement with the TDDE 
results in Fig. 5 of 13] is satisfactory and especially the 
shift of the maximum to higher fields is well reproduced. 
However, in agreement with the discussion in [isj the 
model predicts that the maximum of the ionization rate 
increases monotonically with Z. This is in contrast to 
the numerical findings in [l3| . 



C. Multiphoton ionization 

If many photons are required for ionization, the rela- 
tivistic results may deviate from the non-relativistic ones 
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FIG. 5: Relativistic (blue circles) and non-relativistic (black 
circles) ionization yields obtained for an ion with the nu- 
clear charge Z — 50 exposed to a 20-cycle cos^-shaped laser 
pulse with a peak intensity of 5 x 10^^ W/cm'^ and various 
laser wavelengths. Additionally, the non-relativistic ioniza- 
tion yields (red squares) for an ion with the nuclear charge 
Z' = 50.88 are shown whose non-relativistic ionization po- 
tential is equal to the relativistic ionization potential of the 
ion with Z — 50. The non-relativistic A'^-photon thresholds 
{N = 2 — 4) for Z and Z' are indicated by black and red 
vertical dashed lines, respectively. 



by an order of magnitude or more. This is demonstrated 
in Fig. [5] where relativistic and non-relativistic ioniza- 
tion yields are compared for an ion with the nuclear 
charge Z = 50. The laser wavelength varies in the range 
from 0.05 nm (2-photon ionization) to 0.15 nm (5-photon 
ionization). Whereas the general structure of the wave- 
length dependence is similar for both relativistic and non- 
relativistic treatments, the difference in the positions of 
all pronounced features like peaks and minima may result 
in a substantial discrepancy, if the ion yields are com- 
pared at a single photon frequency. Obviously, the sharp 
steps are due to the A^-photon channel closings and the 
peaks are signatures of resonantly enhanced multipho- 
ton ionization (REMPI) processes. Their positions are 
directly related to excitation energies and the ionization 
potential, respectively, and both are higher in the rela- 
tivistic than in the non-relativistic case. 

In view of the scaling relation (P7| it is interesting to 
compare the TDDE result (for Z — 50) with the non- 
relativistic TDSE result for the scaled nuclear charge 
Z' = 50.88. As can be seen from Fig. [51 the ob- 
tained TDSE results deviate from the relativistic ones 
for Z = 50 only by a few percent. This indicates that 
the change of the ionization potential is the by far dom- 
inating relativistic effect, at least if multiphoton ioniza- 
tion is described within the dipole approximation. Other 
possible effects like the splitting of (resonant) intermedi- 
ate states due to spin-orbit coupling are thus very small, 
but can be quantified more easily, after the scaling re- 
lation has accounted for the main effect. The dominant 
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influence of the sliift of the ionization potential is likely 
to depend on the considered laser intensity and photon 
frequency. It should be reminded that the Keldysh pa- 
rameter 7 — ^J2 IpUj/Fo 36] varies in Fig. [S]in between 
38.17 and 12.72. This is deep in the multiphoton regime, 
in fact even in the perturbative one. 

IV. CONCLUSION 

Single- and multiphoton ionization of highly charged 
atomic ions has been numerically studied by a direct so- 
lution of the time-dependent Dirac equation within the 
dipole approximation. The stationary Dirac equation is 
solved by projecting the radial part onto a _B-spline basis 
and the obtained field-free eigensolutions are used in the 
subsequent time propagation. Results for both length 
and velocity gauges for describing the ion-field interac- 
tion are obtained and compared. The inclusion of the 
negative-energy Dirac states for the description of the rel- 
ativistic dynamics is shown to be important in the case 
of the multiphoton ionization, if the velocity gauge is 
adopted, even if the considered intensities and frequen- 



cies are too low for allowing non-negligible real pair cre- 
ation. If ionization occurs via absorption of a single pho- 
ton or the time propagation is performed in the length 
form, the role of negative energy states is much less sig- 
nificant. 

Comparing solutions of the time-dependent Dirac and 
Schrodinger equations for the same ion and laser pulses, 
the relativistic change of the ionization potential is 
demonstrated to dominate other relativistic effects in the 
case of multiphoton ionization. It is shown that this effect 
is successfully accounted for by a simple scaling relation. 
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